Exact solutions for a rotational flow of generalized 
second grade fluids through a circular cylinder 



Amir MahmoocQ, Saifullah, Qammar Rubab 

Abstract. In this note the velocity field and the associated tangential stress correspond- 
ing to the rotational flows of a generalized second grade fluid within an infinite circular 
cylinder are determined by means of the Laplace and Hankel transforms. At time t = 
the fluid is at rest and the motion is produced by the rotation of the cylinder, around its 
axis, with the angular velocity fit. The velocity field and the adequate shear stress are 
presented under integral and series forms in terms of the generalized G-functions. Fur- 
thermore, they are presented as a sum between the Newtonian solutions and the adequate 
non-Newtonian contributions. The corresponding solutions for the ordinary second grade 
fluid and Newtonian fluid are obtained as particular cases of our solutions for /3 = 1 , 
respectively a = and (5 — 1 . 
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1 Introduction 

The motion of a fluid in a rotating or sliding cylinder is of interest to 
both theoretical and practical domains. It is of very important significance 
to study the mechanism of viscoelastic fluids flow in many industry fields, 
such as oil exploitation, chemical and food industry and bio-engineering 
[1]. Fetecau et al. [2] have considered the general case of helical flow of 
an Oldroyd-B fluid and have determined the velocity fields and the associ- 
ated tangential stresses in forms of series in terms of Bessel functions. Re- 
cently fractional calculus has encountered much success in the description 
of complex dynamics, such as relaxation, oscillation, wave and viscoelastic 
behaviour. Bagley [3], He [4], Tan [5] used fractional calculus to handle 
various problems regarding to flow of the second grade fluid. In this note we 
will study the rotational flow of a generalized second grade fluid within an 
infinite circular cylinder of radius R. The motion is due to the cylinder that 
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at time t = + , begins to rotate around its axis with the angular velocity 

nt. 

Exact analytic solutions of this problem are obtained by using Hankel 
and Laplace transforms and generalized G-functions. Some classical results 
can be obtained as special cases of our solutions. 

2 Governing equations 

The constitutive equation of an incompressible generalized second grade 
fluid is given by [4 - 6] 

T = —pi + /iAi + aiA 2 + a 2 A? , (1) 

where T is the Cauchy stress tensor, —pi denotes the indeterminate spherical 
stress, fi is the coefficient of viscosity, a\ and a 2 are the normal stress moduli 
and Ai and A 2 are the kinematic tensors defined through 

Ai = grad v + (gradv) T , (2) 



A 2 = £>f Ai + Ai (gradv) + (grad v) T Ai . (3) 

In the above relations v is the velocity , the superscript T denotes the 
transpose operator, and Z?f is the Riemann - Liouville fractional derivative 
operator defined by [7] 

where T(-) is the Gamma function. For = 1 the generalized model reduces 

to classical model of second grade fluid because D\f = df jdt. 

Since the fluid is incompressible, it can undergo only isochoric motions and 

hence 

divv = tr Ai = 0. (5) 

If this model is required to be compatible with thermodynamics, then the 
material moduli must meet the following restrictions [8] 

H>0, cci > and a\ + a 2 = . (6) 

In cylindrical coordinates (r,9,z), the rotational flow velocity is given 
by [2, 6] 

v = v(r,t)=u(r,t)ee, (7) 
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where eg is the unit vector in the 9 direction. For such flows the constraint 
of incompressibility is automatically satisfied. 
Introducing (7) into constitutive equation, we find that 

r ( r , t) = (fi + ai D$ )(|- - l) w (r, i), (8) 

where r(r, t) = S r g(r,t) is the shear stress which is different of zero. The 
last equation together with the equations of motion lead to the governing 
equation 

^M = ( v + a Df)(^ + l| : -l) W (r > t), re(0,R), t>0, (9) 

where v = [x/p is the kinematic viscosity, p is the constant density of the 
fluid and a = a±/p. 

3 On the rotational flow through an infinite circu- 
lar cylinder 

Let us consider an incompressible generalized second grade fluid at rest in 
an infinite circular cylinder of radius R. At time zero, the cylinder suddenly 
begins to rotate about its axis with the angular velocity Qt. Owing to the 
shear, the fluid is gradually moved, and its velocity is (7) and governing 
equation is (9). The appropriate initial and boundary conditions are 

w(r,0) = 0; r€[0,R), u(R,t) = Rttt; t>0. (10) 

To solve this problem we shall use as in [6, 9] the Laplace and Hankel 
transforms. 

3.1 Calculation of the velocity field 

Applying the Laplace transform to Eqs. (9) and (10) and using the 
Laplace transform formula for sequential fractional derivatives [7] , we obtain 

where the image function uj(r,q) = J °° uj(r,t)e~ qt dt of oj(r, t) has to satisfy 
the condition 

ZJ(R,q) = ^, (12) 

q being the transform parameter. 
In the following we denote by 

f R - 

u H (ri n ,q)= ru(r,q)Ji(rri n )dr , (13) 
J o 
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the Hankel transform of u(r, q) , where Ji(-) is the Bessel function of first 
kind of order one and r\ n , n = 1, 2, 3, ... are the positive roots of the tran- 
scendental equations J±(Rr) = 0. 

Multiplying now both sides of Eq. (11) by rJ±(rri n ), integrating with re- 
spect to r from to R and taking into account the condition (12) and the 
equality 



R 



doj(r,q) lduj(r,q) uj(r,q) 



dr 2 



+ 



dr 



Ji(rri n )dr 



= Rr ln J 2 (Rn n )uj(R, q) - r ln uj H (r ln , q) , 



we find that 



v + aq 



q + ar\ n qP + vr 



In 



(14) 



(15) 



Now, for a more suitable presentation of the final results, we rewrite Eq. 
(15) in the following equivalent form 



where 



Wff(n„, q) = Lo 1H (ri n , q) + w 2 ff(ri„, q) + oj 3H (ri n , q) , 



£IR 2 

UiH{r ln ,q) = — J2{Rri n ) , 



q"r ln 



nR 2 J 2 (Rr ln ) (\ 



1 



vr 



In 



and 



um{rin,q) = aQR 2 r ln J 2 (Rr lr 



1 



q q + vrf n 



J3-1 



q + vrf n 



q + oir\ n qP + vr 



In 



(16) 

(17) 
(18) 

(19) 



Using the formula 



rR 

I r 2 Ji(rri n )dr = J 2 (-Rri n ) , 

JO r ln 



(20) 



we get that inverse Hankel transform of the function wi#(ri n , q) is the func- 
tion 



(21) 



The inverse Hankel transforms of the functions uikH( r in,q) , k = 2,3 are 
the functions 



UkH{r,q) 



E 



-UkH{rin,q) 



(22) 
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Introducing Eqs. (21) and (22) into Eq. (16) we find that the Laplace 
transform ZJ(r, q) has the form 



Or 2ft ^ Ji(rr ln ) fl 1 * 



a-i 



+ MV ^) 1 2 T £ T . (23) 

To obtain the velocity field to(r, t) = L _1 {uJ(r, q)} we will apply the discrete 
inverse Laplace transform method [6, 7, 9]. For this we use the expansion 



q + arf n qP + vr\ n (q 1 ^ + ar 2 ln ) + vr\ n q-P 

- B-O'^-* Zk+T ' <*> 

fc=o U 1 "/ 3 + arf\ 

Introducing (24) into (23), applying the discrete inverse Laplace transform 
and using the following properties 

L- 1 {F 1 (q)F 2 (q)} = (h * f 2 )(t) = f h{t - s)f 2 (s)ds , 

Jo 

(25) 

where 

fk(t)=L- 1 {F k (q)}, k=l,2, 
L ~ 1 { ( q a 9 l d )c } = G ^Ad,t), Re(ac-b)>0, (26) 

and [10] 

aAc( ' j " ^ r(c)r(i + 1) r[(c + j> - 6] ' (27j 

are the generalized G-functions, we find for c^(r, t) the expression 



OO j i \ oo 

W (r, t) = Ww (r, t) + 2a £ ff £ ( " 



In 



(28) 



x ^ exp[-urj n (t - s)]G 1 - /3i - l3k - ltk+1 ^- ar 2 n ,s^jds , 
where [2, Eq. (4.5)] 

"„(r, = mt - ^ £ r /^ ( ^) [1 - M (29) 

is the similar solution for Newtonian fluids, performing same motion. 
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3.2 Calculation of the shear stress 

Applying the Laplace transform to Eqs. (8) we find that 



d 1, 



r(r, q) = (fi + ai q p ){— - -)v{r, q) . 



(30) 



The image function uj(r,q) can be obtained using Eqs. (27)-(29) and the 
formula 



L {f(^T)} = ^' a> " 1 - 



(31) 



Consequently, applying the Laplace transform to Eq. (28), differentiating 
the result with respect to r and using the identity 

rJ[(rr ln ) - Ji{rr ln ) = -rr ln J 2 {rr ln ) , (32) 

we find that 

doJ uj 2VL ^ J2{rri n ) f I 1 



uuj lu \ <J2\> I In) L 1 

dr r v ^ x r\ n J 2 (Rr ln )\q q + vr\ n 



y r^mn) y 



— v r 



In 



-arj n ) T(k + j + l) 



T(k + l)r(j + 1) 



q + vr 1 Xn q k+(l-/3)(j+l)+l ■ 
Introducing (33) into (30) we get 



r(r,q) = 2 P nY, 



J2(rr ln ) 



(33) 



7/3-1 



^ rf n J 2 (Rr ln ) \q q + vr\ n ) ^ J 2 (i?n n ) 9 + i/r? n 



-2aft^ 



rf n J 2 (rri r 



1 J ^m) £o 



- -w\ n -ar 2 ln ) T(k + j + l) 



T(k + l)T(j + 1) 



A 4 



^air ln 



? + [ g fc+(l-/3)(j+l)+l 9 fc+3+(l- j 9)j-2/3 



g fc+3+(l-/?)j-2/3 



.(34) 
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Applying the inverse Laplace transform to Eq. (34), we find that the shear 
stress r(r, t) has the form 

r(r, t) = r N (r, t) + 2a,Q G, M1 (-< *)" 



.^f i^K) f (-<)'(-^)' r(t+J, + " ;; 
£3 *<*■.) r ( * + i)ro + i, 



X 



<■* „ r Ma *+(i-/9)(j+i) 



/ exp[-urf n (t - s)} 
Jo 



T[k+(1-P)(j + !) + !] 



I/Qir? s fc+2+(l-/9)j"-2/3 ^ °° 2 t / \ 



T[k + 3+(l-P)j-2p]j - ^ J 2 (/2n n ) 



» ( - " ( " « r?„) r ( fc + i + 1) aitfc+2+ (i- j 9) J --2/3 

X J^ r(fc + i)r(j + i) r[fc + 3 + (i-/?)j-2/3] (35) 

where [2, Eq. (5.3) for a = 0] 

r w (r,t) = 2 P nf; [1 - exp(-vrl n t)] , (36) 

is the shear stress corresponding to a Newtonian fluid performing the same 
motion. 



4 Limiting cases 

Making (3 = 1 into Eqs. (28), we obtain the velocity field 

n — 1 k — 

x J exp[-vr\ n {t - s)]G _fc-i,fc+i ^ - arf n , s^j ds , (37) 
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corresponding to an ordinary second grade fluid, performing the same mo- 
tion. Similarly, from (35), we obtain the shear stress 



r(r, t) = r N (r, t) + 2a,Q ^rld, o,i(-^„, *)" 



T(k + l)T(j + l) 



I exp\-vr\ n {t - s)](ji + va x r\ r 
Jo 



T(k + 1 



-ds— 



2aa nf r ^2K) A (-v r\ n ) k {-arljT(k + j + 1) t fc 

1 J2(flrm) fe ^ r(fc + i)r(j + i) r(* + i) 



,(38) 



corresponding to an ordinary second grade fluid, performing the same mo- 
tion. 

The above relations can be simplified if we use the following relations: 



Go_fc_i lfc +i(-ar? n ,s) = f7^—^Yl 



(-arjjrjk + j + l) 
T(k + l)^ Q T{k + l)T(j + l) 



T(k + 1 



-(1 + arl n )- (k+l \ (39) 



°° 1 °° 1 / vr 2 s \ k 

g(-"L)'Go,- l - 1 , +1 (-L,») = TT^Esi- TTi^j = 



1 + arl 



-exp\ — 



In 



vr\ n s 



1 + arl 



In 



and 



G lfi ,i(-vrl n ,t) = exp( -vr\ n t 



(40) 



(41) 



As a result, we find the velocity field and the adequate shear stress under 
simplified forms 



(jj(r, t) = rQt 



2Q Jl( rr lr 



1 — exp 



vr 



In 



1 + arl 



in 



(42) 
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and 



r(r,t) =2pftJ2 



1 + ar\ n 



exp 



vr 



In 



1 



or 



In 



(43) 



which are identical to Eqs. (5.1) and (5.3) from [2]. 

If in Eqs. (42) and (43), we make a = , then the corresponding solu- 
tions of the Newtonian fluids are recovered. 



5 Conclusion 

In this note, the velocity field and the adequate shear stress correspond- 
ing to the rotational flow induced by an infinite circular cylinder in an incom- 
pressible generalized second grade fluid, have been determined using Hankel 
and Laplace transforms. The motion is produced by the circular cylinder 
that at the initial moment begins to rotate around its axis with a constant 
angular acceleration. The solutions that have been obtained, written under 
integral and series forms in terms of generalized G-function, satisfy all im- 
posed initial and boundary conditions. Furthermore, they are presented as 
a sum between the Newtonian solutions and the adequate non-Newtonian 
contributions. In the special case when j3 = 1 , or (3 = 1 and a = , the 
corresponding solutions for ordinary second grade fluid and Newtonian fluid, 
respectively, performing same motion, are obtained. 
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